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•^ ' Abstract. The resources in a cell are finite, which implies that the various 

components of the cell must compete for resources. One such resource is the ribosomes 
used during translation to create proteins. Motivated by this example, we explore 
this competition by connecting two totally asymmetric simple exclusion processes 
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on the effects on the density and current of having different entry and exit rates. 
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^^ • 1. Introduction 

>' 
^ ! Non-equihbrium systems are ubiquitous in nature. With no comprehensive framework 

\0 ! for such systems in general, the understanding of non-equilibrium statistical mechanics 

is recognized as one of the major challenges [ll. To make progress toward finding such a 

CO . 

Q , framework, it is reasonable to study simplified models, in order to gain some insight into 

^ ! this type of complex systems. One such model is totally asymmetric simple exclusion 

process (TASEP). On the one hand, this model is simple enough to be amenable to 
analytic methods, so that many exact results are known. At the same time, it is 
applicable to a wide range of biological and physical systems, e.g. protein production 
[21 El la E], traffic flow P [7], and surface growth [HI E]. 

The simplest version of the TASEP consists of a one-dimensional lattice with 
particles moving unidirectionally from one site to the next. Particles may move 
only if the adjacent site is empty. Two types of boundary conditions are typically 
studied, periodic and open. With periodic boundary conditions, the stationary 
distribution is trivial [lO], though its dynamics differ from that of ordinary diffusion 
[TT| [T2l \T3\ [HI [151 HSl (13 [H]- ^oi open boundary conditions, three distinct phases 
emerge that depend on the entry and exit rates p^ - a low density (LD) phase with 
the lattice less than half filled, a high density (HD) phase with more than half of the 
lattice filled, and a maximal current (MC) phase where the current of particles through 
the lattice is a maximum. If the entry and exit rates are the same, then a shock forms 
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between a LD and HD region that performs a random walk on the lattice. Because 
of the presence of a shock, it is often referred to as the shock phase (SP). The exact 
solution of the steady-state distribution is non-trivial and was found only two decades 
ago [201 |2ll [22]. Not surprisingly, its dynamics is more complex [23| IM] |25| [26l[27[ 128] . 
For a recent review on these aspects of the TASEP, as well as its applications to other 
processes of biological transport, see [29] . 

The TASEP with open boundary conditions has been used to study the production 
of proteins during translation in a cell [21 [3l H]. In this process, ribosomes attach at 
one end of the messenger RNA (mRNA) strand and move unidirectionally to the other 
end. At the other end, the ribosome detaches from the strand and can be used again 
by either the same mRNA or another one. To build more realistic models for protein 
synthesis, modifications to the simplest version of the TASEP have been introduced, such 
as having large particles [3 |30], inhomogeneous hopping rates [3ll [32], and ribosome 
"recychng" [331 (Ml 135]. In this paper, we expand on our previous work on competition 
between multiple TASEPs|ii6j, modeling the simultaneous translation of multiple genes 
in a cell with a limited number of ribosomes. Unlike earlier studies, we consider another 
important aspect of synthesis of proteins in a cell, i.e., the presence of various regulatory 
mechanisms which control the rates of ribosome binding to different proteins. Thus, 
we study TASEP 's with different entry and exit rates. Though we are not aware of 
any similar mechanism for termination, we consider different exit rates also, simply as 
part of a systematic investigation. With such a large parameter space to explore, we 
restrict ourselves to only two TASEPs here, in search for novel and (possibly) universal 
properties that could be applicable for mRNA competition in a real cell. 

This paper is organized as follows: In the next section, we define our model. In 
section [3l we present our simulation results. We give some theoretical considerations in 
section |H Finally, we give a summary and outlook in section [5] 

2. Model specifications 

In our previous study [36] , we model the competition between mRNAs by couphng two 
or more open TASEPs to a finite pool of Np particles and let the entry rates depend 
on this Np. Particles exiting each TASEP join this pool and are "recycled" for entry 
into any of the other TASEPs. Thus, the total number of particles Ntot is conserved. 
While on any lattice, the particles move uni-directionally from one side to the other as 
in the ordinary TASEP. All internal hopping rates are set to unity. Our current model 
(shown in figure [1]) differs from [36]: Here, we relax the constraint that the intrinsic 
(i.e., limiting) entry rates of the TASEPs are identical. Thus, we define ai^2 as the 
intrinsic rate for our two- TASEP system, applicable when the supply of particles is very 
large. For simplicity, let us assume the crossover function (/) to be the same, so that 
the effective entry rates are given by 

Oieff,l = OilfiNp) (1) 

ae//,2 = a2f{Np) (2) 
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Figure 1. Our current model of connecting two TASEPs to a finite pool of particles. 
The large arrows indicate how particles enter and leave the pool. 



As in E 



we will use 



f{N,) = tanh 



N* 



(3) 



(where N* is a scaling parameter), so that / (0) = and / — )■ 1 as A^p — )■ oo. Clearly, it 
is reasonable to use the labels "faster" / "slower" TASEP for the one with larger/smaller 
a. We also consider different exit rates /3i,2, even though we are not aware of biological 
systems which exhibit such differences. 

In our Monte Carlo simulations, we first consider the case of two TASEPs of lengths 
Li and L2 connected to a single pool of particles. To represent the pool, we have a 
"virtual" site, with unlimited occupation (so that we have Li + L2 + 1 sites in total). 
Since this site is connected to both TASEPs, there are actually Li + L2 + 2 "bonds" 
connecting the sites. The simulations are performed as follows. In an update attempt, 
we randomly choose one bond and update the contents of the sites according to the usual 
rules: A hole-particle pair within a TASEP is left unchanged, while a particle-hole pair 
is always changed to a hole-particle pair. If a pool- TASEP bond is chosen and the entry 
site is empty, then a particle is moved in it with probability ae//,i or ae//,2- Finally, for 
the TASEP-pool bond, a particle in the last site is moved into the pool with probability 
/9i,2- One Monte Carlo step (MCS) is defined as Li + L2 + 2 attempts. 

Starting with Nfot particles in the pool (none on the TASEPs), we allow the system 
to reach steady-state, which typically takes 100k MCS. For the next IM MCS, we record 
the density profile (p (x)) for each TASEP at every 100 MCS. From these, we compute 
the overall densities (p), for a total of 10k data points. We also measure the average 
currents (J), by measuring (for example) the total number of particles which exit each 
TASEP over the run and dividing that by 10^. As in the earlier study, we are interested 
in how these quantities are affected by varying Ntot- The profiles obviously contain much 
more detailed information. Thus, in this first stage, we will mostly report the behavior 
of the four functions pi^2 (Ntot) and Ji_2 (Ntot)- 
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Our model has a total of eight parameters: Li, L2, ai, 02, /^i, (^2, Ntot, and A^*.To 
reduce the number of parameters, we fix A^* = 1000. A^* controls the strength of the 
feedback effect for both TASEPs; however, we are focusing on the effects of having 
different entry and exit rates, so we will not explore the effects of A^* in this study. 
Since we have different a's and /3's, each TASEP can be in a different phase (LD, HD, 
MC, or SP) when the pool size becomes large. Thus, 16 different combinations are 
possible. From our experience [M| [36] , the most interesting phenomena occur in the 
combination HD-HD, the results of which will be presented next. 



3. Simulation Results 



3.1. HD-HD 



From the earlier study \3^, the overall density of a constrained HD-TASEP displays 
three regimes, as Ntot is increased: an LD dominated one, a "crossover regime", and one 
controlled by HD. Respectively, these are characterized by aeff (Np) < (3 , «£// (Np) = f3, 
and aeff (Np) > /3. In the crossover regime, Np remains fixed, while all changes in Ntot 
are absorbed by the lattice. Thus, p increases linearly, from the LD value of /3 to the 
HD value of 1 — /3. The threshold values of Ntot are given by aeff {Ntot — PL) = f3 and 
aeff {Ntot ~ {^ ^ f^) L) = /3. These characteristics are again present when two TASEPs 
compete for the pool. The novel features here are the following. If the two TASEPs 
make their crossovers at entirely different points, then all changes in Ntot are absorbed by 
whichever is in the crossover regime, so that activity in both the pool and its competitor 
are completely interrupted. In Figure [2], we illustrate this phenomenon with the case of 
Li = 1000, L2 = 1000, ai = 0.8, Pi = 0.2, ^2 = 0.6, and P2 = 0.4. Note first that the 
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Figure 2. Two TASEPs of equal lengths and different rates with a's and /3's in the 
HD phase. The open circles and squares are the results from the domain wall theory 
presented in section [ 



two TASEPs fill at different rates at low Ntot- This difference is a simple consequence 



of ae//,i — aiNp/N* > a. 



eff,2 



a2Np/N*. Next, from Ntot ~ 600 to ~ 1200, the 
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faster TASEP makes its crossover while the numbers in the pool and the slower TASEP 
remains constant. Thereafter, the slower TASEP continues on its LD regime and, lastly, 
makes its crossover in, approximately, the interval [2000,2200]. We emphasize that, in 
the respective crossover regimes, pi G [/3i, 1 — (3i] and p2 E [(32,1 — (32]- 

To understand this effect, we examine the density profile. Even after the faster 
TASEP reaches the HD state, its entry rate continues to increase. This increase results 
in the decay of the tail near the entrance to changing as Ntot increases, similar to 
changing a (with fixed /3) in the unconstrained, ordinary TASEP |20l EH [22]. As the 
slower TASEP moves through a crossover regime, the tail in the profile of the faster 
TASEP does not change. During each crossover from LD to HD, the average number 
of particles in the pool remains constant. Since each Oe// depends on Np, the ae//'s 
also remain constant as Ntot increases. The extra particles from the increase in Ntot 
are added to the TASEP crossing the phase boundary between the LD and HD phases, 
resulting in the formation of a localized shock. A similar phenomenon is found in a 
single constrained TASEP [3^ and multiple TASEPs with the same a and (3 [35] . 

For both TASEPs to be in the crossover regime simultaneously, each Oe// must 
reach /3 at the same Ntot value. This condition is achieved when ai/f3i = a2//32- The 
slower TASEP 's overall density increases linearly with Ntot in the crossover regime, 
but the faster TASEP's density does not. Two examples are shown in figure |3] for 



ai = 0.8, I3i = 0.2, a2 = 0.6, /32 = 0.15 and (b) 
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LO, (3i = 0.4, 



^2 



L2 = 1000, 

0.5, /32 = 0.2. The ratio of ai/a2 controls the formation of a plateau region for 
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Figure 3. Overall densities and currents of two TASEPs of equal lengths entering the 



crossover regime simultaneously with (a) ai/a2 = 1.33 and[(b)]ai/a2 = 2 



the faster TASEP with a density of p = 0.5. In our previous study [36], similar plateau 
regions form when the lengths of the TASEPs differed. 

Another way to visualize the difference of having both TASEPs in the crossover 
regime is by looking at the probability P{Ni,N2) to find the A^i particles on the first 
TASEP and A^2 particles on the second. When ai/(3i ^ 0.2/ h and in the crossover 
regime, the distribution is sharply peaked about the average Ni and A^2; otherwise, it is 
spread across a range of particle occupation pairs whose sum is constant. These cases 
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are shown in figurelUfor (a) Li = L2 = 1000, ai = 0.8, /3i = 0.2, 02 = 0.6, and /32 = 0.15 
with Ntot = 1250 and Li = L2 = 1000, ai = 0.8, /3i = 0.2, 02 = 0.6, and P2 = 0.40 with 
(b) Ntot = 900 and (c) Ntot = 2100. The increase of the spread of the distribution when 




Figure 4. Distributions of particle occupation for two TASEPs when (a) both TASEPs 
are in the crossover regime, (b) one TASEP in the crossover regime and the other in a 
LD state, and (c) one TASEP in the crossover regime and the other in a HD state. 



both TASEPs are in the crossover regime comes from the additional degree of freedom 
that the second TASEP provides in keeping the average number of particles in the pool 
constant [32] • It is important to note that the ranges of A^ values are the same for both 
TASEPs and are governed by the exit rate of the faster TASEP, A"! 2 G [f3L, (1 — P)L]. 
To further investigate this crossover regime, we turn to the density profile. Here, 
we find that the confinement of the shock between the LD and HD regions is controlled 
by the ratio of ai/a2. Figure [5] shows the density profiles for the same set of parameters 
shown in figure [3] at f(a)] Ntot = 1250 and f(b)] Ntot = 1400. For the simplest TASEP 





(a) (b) 

Figure 5. Density profiles of two TASEPs of equal lengths entering the crossover 



regime simultaneously for (a) ai/a2 — 1.33 with Ntot — 1250 and |(b)| Q:i/Q:2 = 2 with 
Ntot = 1400. 



in the SP, the shock performs a random walk over the entire lattice, which results in a 
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linear density profile ^U\ [211 E2] • The linearly increasing regions in the profiles in figure 
[5] indicate the allowed portions of the lattice on which each shock performs a random 
walk. The fiat regions (of LD or HD) are areas in which the shock does not travel. 

Since the number of particles in the pool remain relatively constant in the crossover 
regime, excess particles are free to choose which lattice to occupy. Due to the constraint 
of ai/(3i = a2//32, the ai/a2 ratio correlates with the difference between the two shock 
heights (i.e. the difference between the LD and HD regions densities). The faster TASEP 
will always have a smaller shock height, which limits the range of the number of particles 
it can hold, [(3L, (1 — f3)L]. The same particle limit applies the the slower TASEP as 
shown in figure HI But due to larger shock height, the shock is now confined to a smaller 
portion of the lattice than the faster TASEP's shock in order to have the same range of 
A^ values (thereby keeping the pool size relatively constant). By decreasing the shock 
height in the faster TASEP, the range of particles it can hold decreases. Thus, the shock 
becomes confined over a smaller region on the slower TASEP as seen in figure |5l This 
effect was not seen in our previous study [36] since it is a result of having different entry 
and exit rates. 

When both TASEPs enter the crossover regime at the same Ntot and the lengths 
are unequal, we see a trend in the overall density similar to the case of equal rates in 
[36] . shown in figure |6] for Li = 1000, L2 = 100, ai = 0.8, /3i = 0.2, ^2 = 0.6, and 
P2 = 0.15. The smaller TASEP has a density of 0.5 for most of the crossover regime. 
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Figure 6. Two TASEPs of unequal lengths with different a's and /3's entering the 
crossover regime at the same Ntot value. 



quickly rising to this value from the LD state and from this value to the HD state. The 
density profile for this TASEP is linear indicating a delocalized shock. Figure [7] shows 
this derealization for Li = 1000, L2 = 100, ai = 0.8, /3i = 0.2, 02 = 0.6, P2 = 0.15, 
and Ntot = 800. The larger TASEP has a localized shock during the crossover regime 
as seen in the density profile in figure [3 Even when the rates are reversed, the smaller 
TASEP has a delocalized shock. We can conclude that, as long as the size of the smaller 
TASEP is less than the intrinsic width of the shock localization, the smaller TASEP 
will have a delocalized shock. 
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Figure 7. Density profiles for two TASEPs witli different lengths. 



3.2. HD-SP, HD-LD, and HD-MC 

The combination of having a and (3 on one TASEP in a HD phase with a and (3 on the 
other TASEP in another phase produces an effect on the density and current similar 
to having different ratios of a/P for each TASEP. Initially, both TASEPs are in the 
LD state when Ntot is small. As we increase the number of particles in the system, 
the HD TASEP begins to crossover from the LD state to a HD one, while the other 
TASEP 's density and current remain constant during this regime. After the HD TASEP 
enters its HD state, the other TASEP's density and current continue to increase until 



it reaches its final state. Examples of this effect are shown in figures 8(a), 8(b), and 
8(c)] for Li = L2 = 1000, ai = 0.7, /3i = 0.3 and ^2 = /32 = 0.3, a2 = I - P2 = 0.3, 
«2 = /32 = 0.7, respectively. We see no new phenomena when we have different lengths 
for the TASEPs. Further, as figure M shows, an appropriately generalized domain wall 



(GDW) theory, which is presented in section 14. 4[ is quite adequate in predicting the 
behavior of the overall density (and therefore the current as well). 



Finally, we present the data for the HD-MC combination (figure 8(c) ). The various 
regimes here are easy to understand qualitatively. Since ai = 02, the initial rise of the 
densities are the same. Thereafter, if there were no competition, the behavior of the 
second TASEP would rise smoothly until a2,e// reaches 0.5. But this rise is interrupted 
by the first TASEP traversing the crossover regime, i.e., in the second region. In the 
third region, it continues its increase and, in the last section, it remains in the MC 
phase. As the first TASEP has essentially dropped out of the competition, there is no 
"kink" in the transition between these regimes, of course (as better displayed by the 
current around A^tot ~ 2000). 



3.3. Other phase combinations 

When the a's and /3's are such that neither TASEP will enter the HD phase, we find the 
density and current behaving in a manner similar to a single, constrained TASEP [34] 
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Figure 8. Two TASEPs with one in the HD state and the other in the (a) SP, Kb 



LD, (c) MC state. 



as the pool size is increased. While the total number of particles needed to saturate the 
system is larger than the number needed for a single TASEP, we find no new features 
emerging in the overall density and current as a function of Nfot, even for different 



lengths. Some typical results are shown in the figures found in [Appendix A 



4. Theoretical considerations 

While we presented some qualitative analysis in the previous section, we now supplement 
those results with a more quantitative analysis for the various phase regimes. 



4.1. LD state 

Regardless of the entry and exit rates, both TASEPs are in a LD state when Ntot is small 
when compared to the smallest lattice length. From the ordinary TASEP [211 EH 138] . 
we know that the overall density is given by p = a; and for a single TASEP with finite 
resources [51] it is equal to the average effective entry rate, p = aeff- Extending these 
results for two TASEPs with unequal entry and exit rates, we have 



Pi = «e//,l = otif{Ntot - PiLi - P2L2) 
P2 = ae//,2 = a2f{Ntot - PiLi - P2L2) 



(4) 
(5) 
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The tte/z's depend on both pi and P2; therefore, a self-consistent solution is found using 
these two equations. For the / chosen in this paper, the solution is found numerically. 
Once one of the TASEPs has left the LD state, we must modify our equations. 



4.2. MC state 

If one of the TASEPs enters the MC state, an increase in a^ff no longer has an effect on 
the density or current. The transition occurs when its Oe// = 1/2. As with the ordinary 
TASEP |20l [211 122], the density p = 1/2 and the current J = 1/4. The Nm at which 
this TASEP reaches is final density (assuming pi is entering the MC state) is 

r' (7^) + ^ + P2L. (6) 



Ni 



tot 



2ai 



where p2 is the density of the second TASEP. This density depends on its state, 

P2=S -2 MC (7) 

I-/32 HD 



For the parameters shown in figures A 1(e) and A 1(b) (where the other TASEP is in 
the LD state), the MC state is reached at N^ot ~ 1600. In figure 8(c), the MC state 
is reached at Ntot ~ 2100 with the other TASEP in the HD state. Both TASEPs are 



approaching the MC state as Ntot increases in figure A 1 (d) , where the first one reaches 
its final density at Ntot ~ 1600 and the second one at Ntot ~ 2200. These values agree 
with the values predicted by equations (Q and ([7]). 



4-3. HD crossover 

The HD TASEP enters the crossover regime as Ntot increases when «£// = P = p and 
leaves when p = 1 — /3 \3^. Taking the HD TASEP to be pi, the beginning Ntot,i value 
of the crossover is 

where the value of p2 depends on the state of the second TASEP. Two possibilities exist: 
the second TASEP is in either the LD state or HD state. Then p2 is given by 

a\ 

P2 



N 



tot,l 



r 



+ (3iLi + P2L2 



(8) 



P2 



a2^ LD 



1 



HD 



(9) 



For the parameters shown in figures 8(a) and 8(b) with the second TASEP in the LD 



state, equations (|8]) and ([9]) give a value of Ntot,i — 887, which agrees with the simulation 



results. Using the parameters for figure 8(c) , we obtain Ntot.i — 1058, which also agrees 
with the data shown. We have similar agreement for figure [2] with Ntot,i — 705 and 
Ntot,i ^ 1905 for each TASEP. 

When leaving the crossover regime, the Ntot value is given by 



N 



tot,2 



r ('^ 



«! 



;i-/3i)Li + p2^2 



(10) 
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where p2 is given above. Equations ( ITU]) and (^ give Ntot,2 — 1286 and Ntot,2 — 1458 for 



the parameters shown in figures 8(a) and 8(c), respectively. For the parameters shown 



in figure O equations ( TTOl) and ([9]) result in Ntot,2 — 1205 and Ntot,2 — 2204 for each 
TASEP. All these values agree with the simulation results. 

Determining the density of each TASEP during the crossover regime is simple when 
only one of the TASEPs is in this regime. The density of the one in the crossover regime 
rises linearly with Ntot, similar to a single constrained TASEP [31], while the other 
density remains constant. Taking pi to be crossing over, we have 



.-1 / A 



piL, = Ntot-r'{^j +P2L2 (11) 

where p2 is in either the LD or HD state as before. However, this simple approach does 
not work if both TASEPs enter the crossover regime at the same time. To understand 
how the density varies with Ntot in that situation and with the SP case, we turn to a 
domain wall approach. 

4.4- Domain wall theory 

The phenomenological domain wall theory has been successfully applied to the 
unconstrained TASEP [391 SOI E] as well as ones with finite resources [351 EH] to 
understand the steady-state results. The theory assumes the presence of a sharp domain 
wall, or shock, separating a low density region near the entrance of a TASEP and a high 
density region near the exit. The shock's movement on the lattice depends on the 
currents of the particles (holes) entering from the entrance (exit) and the wall height 



For an ordinary TASEP with no feedback mechanism, the domain wall moves to 
the left and to the right with fixed rates that depend on a and /3 pTj. We generalize 
this result to include the feedback effect of Oe//; thus, the hopping rates become site 
dependent. While we cannot use the generalized domain wall (GDW) theory when either 
TASEP is in the MC state, we apply it to all other cases here. Due to the connection 
between the shock positions /ci, /c2, and Np, 

Np = Ntot-{^-f3i){L^-h)-aJ{Np)h-{l-f32){L2-h)-a2f{Np)h{12) 

the function f{Np) can be rewritten as f{ki,k2). Then / is found using this self- 
consistent equation. As the values of ki and ^2 (and subsequently /) change, the 
current of incoming particles ae//(l — cte//) and domain wall height 1 — f3 — aefj for 
each TASEP will also change due to their dependence on /. This fluctuation in / will 
lead to domain wall hopping rates that are site dependent [351 ES] . 

With two TASEPs connected to a single flnite pool of particles, the probability P 
of finding a set of domain wall positions {ki,k2} at steady-state is given by 

= Dl_,^,^P{k, - 1, k2) + /^fc- +i,fe,P(/^i + 1, A:2) + El,,^_,Pik^, k2 - 1) 

+E^,M+iP(^^^ A;2 + 1) - {D^,, + /^fc- ,,, + El,^ + E,- _,J P{k„ k2) (13) 
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where 



D 



kiM 



D 
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E 



c^ifjh, k2){l - aifih, fca)) 
1 -/3i -aif{ki,k2) 
/3i(l-/3i) 
l-/3i -aif{ki,k2) 
"2/(^1, ^2)(1 - a2f{ki, h)) 
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1-/32 
/3i(l- 



- a2fiki,k2) 

P2) 



kiM 



(14) 
(15) 
(16) 

1-/32 -oi2f{ki,k2) 

along with appropriate reflecting boundary conditions. We lose the detailed balance 
that was previously exploited to find an analytical solution |36]. While it is possible to 
find the P{ki, ^2) analytically, it is not very practical. This system of (Li + 1)(L2 + 1) 
equations ((Li + 1)(L2 + 1) — 1 which are linearly independent) becomes time-consuming 
to solve even for numerically finding the eigenvector corresponding to the zero eigenvalue. 
Instead, we build the probability distribution through Monte Carlo simulations of a 
random walker on a two-dimensional lattice with the hopping rates D^, D~ , E^ , and 
E~ . These simulations give us the P(fci, k2) we need to calculate the density profile and 
overall density. The profile for each TASEP is given by [36] 

L2 



Pii.x) = 5Z 



^ (1 - /3l)P(fcl, k2) + Y, «l/(^l' ^2)P(fcl, k2] 



p2[X) 



k2=0 Lfci=0 
Li 

E 

fci=0 



X 

E( 

,fc2=0 



x+1 
L2 



(5,)P{ku k2) + Y, «i/(^i' ^2)P(A;i, k2) 



x+l 



(19) 



The overall density is given by pi = J2x=iPii^)- "^^^ GDW theory results agree with 



the simulation results as shown in figures [21 [3l O 8(a), 8(b), A 1(a), A 1(c), A 1(f) for the 



overall density, and figures [5l [3 for the density profile. 

The domain wall picture helps explain the difference between the results in figure 



[6] and 3(b) In figure [6], the delocalization of the shock over a range of N^ot is due 



to the domain wall refiecting at the boundaries on the smaller TASEP. In figure 3(b) 



the difference in hopping rates allow the shock in the faster TASEP (larger rates) to 
move about the entire lattice more easily than the one on the slower TASEP (smaller 
rates). The shock in the slower TASEP will be less likely to move away from its average 
position, leading to shock localization. Also, the domain wall height, which appears in 
the denominator of the hopping rates, plays a significant role. If the wall height is too 
large, then the difference between the rates for each TASEP decreases. The smaller 
difference allows the shock to wander over a large portion of the slower TASEP, as seen 



in figure 5(a) Thus, shock localization can be induced by either different lengths or 
different rates. 



Finally, associated with figure 8(c) (HD-MC), we have no GDW theory to provide 



a good theoretical prediction, as the second TASEP ends in a state with no domain 
walls (MC). Since the general aspects of this competition is qualitatively understood. 
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designing a more sophisticated and quantitative theory seems unnecessary. 

5. Summary and Outlook 

In this paper, we explored how competition for particles between two TASEPs affect the 
overall density, density profile, and current. Through simulation results and theoretical 
considerations, we have shown that new effects arise from having different entry and exit 
rates on the TASEPs. One of these effects is the localization of a shock on the lattice 
due to the difference in entry and exit rates. The appropriately generalized domain wall 
theory captured the shock localization phenomenon and reproduced the overall density 
and density profiles. However, more work still needs to be done if we want to make a 
connection to the translation process in a cell. 

While our study has focused on only two TASEPs, more should be added. Recalling 
our motivation of protein synthesis, many mRNA's compete for the same pool of 
ribosomes. The parameter space to explore increases with each additional TASEP, 
which could lead to new phenomena occurring. Similarly, the dimension of the random 
walk set forth in the GDW theory increases for each new TASEP that is added. A 
systematic study of multiple TASEPs would be useful. 

Beyond multiple TASEPs, other additions to the model should be made in order 
to better model the translation process during protein synthesis [21 [Sj SJ |5l |12]. First, 
the ribosome does not move to the next codon at the same rate for all codons, and 
the rate may depend on the concentration of amino acid transfer-RNAs (aa-tRNA) in 
the cell [13]. Thus, TASEPs with inhomogeneous, mRNA-sequence dependent, hopping 
rates must be taken into account [11I32]. Now that these rates depend on the aa-tRNA 
concentrations, it is reasonable to consider the competition for finite aa-tRNA resources. 
Notably, such an ambitious undertaking has been carried out recently [^ H5] , although 
the behavior in a real cell, with thousands of copies of thousands of different genes, 
will remain difficult for simulation studies in the conceivable future. Second, ribosomes 
cover more than one codon, typically 12 [16]. Therefore, the size of the particles should 
be larger as well [21 |3l IH El ^0\. By combining these individual elements, we hope to 
gain a better understanding of the translation process during protein synthesis, as well 
as non-equilibrium systems in general. 

After completing this work, we became aware of a similar study by P. Greulich, 
et. al. [1?]. The main differences between our efforts are the following. 1) We explore 
the density profile in our Monte Carlo simulations and theoretical approaches. 2) We 
distinguish between systems in the SP with localized shocks and those with delocalized 
ones. 3) We explain our results from a domain wall perspective for both the overall 
density and density profile, instead of using a mean-field approach as in [17] . 



Competition for finite resources 14 

Acknowledgments 

We would like to thank Jiajia Dong and Beate Schmittmann for insightful discussions, 
Irina Mazilu and Tom Williams for a critical reading of the manuscript, and Martin 
Evans for calling our attention to ref. j47]. This work was funded in part by the U.S. 
National Science Foundation through Grant No. DMR-1005417, and Washington and 
Lee University through the Lenfest Grant. 



Competition for finite resources 

Appendix A. Results for cases without an HD phase 



15 



Results in figurelAllare shown for Li = L2 = 1000, A^* = 1000, and (a) ai = 1— /3i = 0.3, 
a2 = 1- P2 = 0.4; m\a, = 1- P, = 0.3, aa = /^s = 0.7; M\ 



ai 



l-/3i 



0.3, 



aa = /32 = 0.3;[(d)]ai = /3i = 0.8, as = /Sa = 0.6;[(e)]ai = /3i = 0.7, aa = /92 = 0.3; [(f) 
ai = /3i = 0.3, aa = /32 = 0.4. 
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Figure Al. Overall density and current for various phase combinations. 
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